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Generalized parton distributions (GPDs) are studied at the hadronic (nonperturbative) 
scale within different assumptions based on a relativistic constituent quark model. In 
particular, by means of a meson-cloud model we investigate the role of nonperturbative 
antiquark degrees of freedom and the valence quark contribution. A QCD evolution of the 
obtained GPDs is used to add perturbative effects and to investigate the GPDs' sensitivity 
to the nonperturbative ingredients of the calculation at larger (experimental) scale. 
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1. Modeling GPDs with double distributions 
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^ ■ We will concentrate our attention on the chiral even (helicity conserving) distribution 
H q (x, £, Q 2 , t) for partons of flavor q at the hadronic scale where the models we are going 
to discuss are assumed to be valid to evaluate the twist-two amplitude. The invariant 
momentum square is t = A 2 = (P /M — P^) 2 , x is the quark light-cone momentum fraction 
with respect to the average nucleon momentum P M = |(-P M + P" 1 ), and the skewedness 
£ describes the longitudinal change of the nucleon momentum, 2£ = — A + /P + . In the 
following the dependence on the scale Q 2 is understood. 

We also introduce non-singlet (valence) and singlet quark distributions, 

H NS (x, £,*) = E [H«(x, £, t) + H*{-T, £, t)] = H NS (-x, £, t), (1) 
q 

H s (x, l Hq & Z> - HQ (-^ Z> *)] = ~H s {-% £, t), (2) 

respectively. Besides being symmetric or antisymmetric in x, they are also symmetric 
under £ — >• — £ due to the polynomiality property [|2]. The analogous GPD for gluons is 
symmetric in x, i.e. H 9 (x,£,t) = H 9 (—x,^,t), and reduces to the gluon density g(x) in 
the forward limit (x — > x, H 9 (x, 0, 0) = xg(x), x > 0). 

Following ref. [ H] we will assume a factorized t dependence determined by some form 
factor, and parametrize the t-dependent part in terms of double distributions (DDs) in- 
volving a given profile function and the forward parton distribution q(x) derived in some 
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model. At the hadronic scale Q$, where the short range (perturbative) part of the in- 
teraction is negligible and, therefore, the glue and sea are suppressed, the long range 
(confining) part of the interaction produces a proton composed by (three) valence quarks, 
mainly [ 3 r Therefore quark models are suitable to construct the parton distribution at 
the scale Ql [Ej. 

According to the method discussed in two recent papers [ EJ E] along the lines of 
ref- [Ej, we evaluate the valence contribution q bare (x) to the parton distribution within 
relativistic light-front constituent quark models (CQMs) at the scale Qq. This distribution 
automatically fulfills the support condition and satisfies the (particle) baryon number and 
momentum sum rules. Results presented here are obtained for the hypercentral CQM [ 

El. 

2. Parton distributions and the meson-cloud model 

The meson-cloud model introduces the sea by incorporating qq contributions into the 
valence-quark model of the parton distribution discussed in the previous section. The 
basic hypothesis of the meson-cloud model is that the physical nucleon state can be 
expanded (in the infinite momentum frame (IMF) and in the one-meson approximation) 
in a series involving bare nucleons and two-particle, meson-baryon states [Ej. We will 
consider fluctuations of the proton including both pion-nucleon and pion-Delta states. 
The partonic content of the A and the pion will be consistently evaluated within the 
same scheme assuming light-front dynamics and valence contributions only. 

The quark distributions in a physical proton are then given by 

q p (x)=Zq b p are (x) + Sq p (x), (3) 

where Sq p (x) is the contribution coming from the meson-baryon fluctuations, and Z is a 
suitable renormalization constant. 

Matching the sea, valence and gluon distributions within the radiative parton model 
we start evolution with continuous functions all over the range —1 < x < 1 and identify 
the matching scale Ql = 0.27 GeV 2 consistent with QCD evolution equations [ ITUl ITT]. 

3. Results and discussion 

Results are presented in this section according to the model based on DDs. The t 
dependence is dropped from the very beginning and could be reintroduced in the final 
results by an appropriate t-dependent factor. The D-term is omitted as well. The QCD 
evolution was numerically performed to NLO accuracy according to a modified version of 
the code of ref. [H2j (see [H3] for further details and results). 

The singlet quark, non-singlet quark and gluon GPDs obtained in the model have 
been studied as a function of x, £ and Q 2 . With no initial gluons and an input parton 
distribution given by q bare (x) the model already gives a nonvanishing contribution to quark 
GPDs in the ERBL region \x\ < £ at the hadronic scale without introducing discontinuities 
at x = £ and with a weak £ dependence. In particular, the absence of the sea contribution 
gives H s = H NS at x > £. After evolution up to Q 2 = 5 GeV 2 GPDs are almost confined 
into the ERBL region with a significant £ dependence. 
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Figure 1. Singlet and non-singlet quark and gluon GPDs at £ = 0.2 obtained from 
NLO evolution within the double-distribution model using parton distributions of the 
hypercentral CQM with the sea contribution at the initial scale Qq = 0.27 GeV 2 . 
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In fig. [T] results are shown for QCD evolution up to Q 2 = 5 GeV 2 of the singlet quark, 
non-singlet quark and gluon GPDs obtained with the parton distributions of the hyper- 
central CQM and including the sea contribution at the initial scale Qq = 0.27 GeV 2 . The 
distributions are plotted at £ = 0.2 as a function of Q 2 . In fact, the largest effects of 
evolution modify the input GPDs within the first few GeV 2 in the Q 2 evolution, Q 2 = 5 
GeV 2 being a value where GPDs have already reached a stable configuration with respect 
to their asymptotic shape. 

The model already gives a nonvanishing gluon contribution at the hadronic scale. Under 
evolution as the resolution scale increases the distributions are again swept from the 
DGLAP domain to lie almost fully within the ERBL region. This is a consequence of the 
fact that functions with support entirely in the time-like ERBL region are never pushed 
out of the ERBL domain. In fact, the evolution in the ERBL region depends on the 
DGLAP region, whereas the DGLAP evolution is independent of the ERBL region. The 
qualitative result after evolution is similar to the case without the sea [ ITH] in spite of a 
more pronounced £ dependence of H s at the input hadronic scale and a shape of H NS 
sensitive to the input at x — 0. 

The present results focus on the ERBL region as the most interesting one to look at the 
nonperturbative effects surviving evolution. This is suggesting that one has to investigate 
suitable processes under appropriate kinematic conditions to study such effects. 
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